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Abstract 

We derive a relativistic expression for neutrino energy losses caused by the 

direct Urea processes in degenerate baryon matter of neutron stars. We use 

two different ways to calculate the emissivity caused by the reactions to our 

interest. First we perform a standard calculation by Fermi's "golden" rule. 

The second calculation, resulting in the same expression, is performed with 

the aid of polarization functions of the medium. Our result for neutrino 

energy losses strongly differs from previous non-relativistic results. We also 

discuss nonconservation of the baryon vector current in reactions through 

weak charged currents in the medium, when the asymmetry between protons 

and neutrons is considered. The above effects, not discussed in the literature 

before, substantially modify the polarization functions responsible for the 

induced weak charged currents in baryon matter. 
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In numerical simulations, neutrino cooling of the massive core of a neutron star, with 
the standard nuclear composition, is governed by the direct Urea processes n p + e~ + Ue, 
p -\- ^ n + Ue- The important role of these reactions in the rapid cooling of neutron 
stars has been first pointed out by Boguta [|T|. Later Lattimer et al. [0] have suggested a 
simple formula for neutrino energy losses caused by the above processes in degenerate nuclear 
matter under /^-equilibrium, which exhibits a threshold of proton concentration necessary 
for the direct Urea processes to operate^. The above estimates have been made in a non- 
relativistic manner, assuming that the participating particles are free. Some improvement, 
which partially accounts for strong interactions, was made by replacing the baryon masses 
with their effective values. 

Actually, the superthreshold proton fraction in the core of neutron stars appears at large 
densities, when Fermi momenta of participating nucleons are of the order, or larger, than 
their effective masses. The total four-momentum of the final lepton and antineutrino is 
time-like; therefore, in the free relativistic gas, the energy-momentum conservation requires 
a large difference in the effective masses of protons and neutrons M* — M* ~ lO'^MeV that is 
unlikely to appear. Thus, in the relativistic regime, the energy conservation can be fulfilled 
only by taking into account the difference in the potential energies of proton and neutron. A 
further simplification was made by neglecting the proton recoil. This is not justified in the 
relativistic regime, because the momentum of the final lepton is of the order of the proton 
effective mass. Both the proton recoil and the difference in the proton-neutron potential 
energies strongly modify the emissivity of the direct Urea reactions even in the Mean Field 
Approximation 

In the present paper we derive a relativistic expression for neutrino energy losses caused 
by the direct Urea processes by taking into account the above effects. We use two different 



^Prakash et al. Q have found that an admixture of A hyperons opens the direct Urea processes 
in baryon matter for an arbitrary proton concentration. 
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ways to calculate the emissivity caused by the reactions to our interest. First we perform a 
standard calculation by Fermi's "golden" rule. The second calculation, resulting in the same 
expression, is performed with the aid of polarization functions of the medium. In addition 
to allow us to check the result obtained by the previous method, this second method shows 
some important details, which must be taken into account in the next step (when RPA 
corrections are included). 

Relativistic polarization functions of hot baryon matter have been investigated earlier 
by several authors (see e.g. [Q, P], and references therein) having made analytic com- 
putations for the case when the particle-hole excitations correspond to identical baryons. 
The study of interactions through charged weak-currents requieres considering polarization 
functions of another type, in which the virtual particle and the hole correspond to baryons 
of different kinds. As we show below, the difference in potential energies of participating 
particles leads to non-conservation of the vector transition current of baryons in the direct 
Urea reactions. This effect dramatically modifies the corresponding polarization functions 
of the medium and must be taken into account in order to obtain the correct result. 

According to mean-field theory, the medium-modified energy of a relativistic baryon 
of kind B is similar to the energy of a single particle with effective Dirac mass in 
the effective potential Ub of self-consistently generated meson fields^. The effective mass, 
as well as the effective potential, arises due to interactions of the fields. We denote by 



(p) = yP^ + Mff the baryon kinetic energy. The single-particle energy of a baryon with 
momentum p is Eb (p) = Eb (p) + Ub, so that the individual Fermi distributions are of the 
form 



where /is is the chemical potential of 5- kind baryons. 

The low-energy Lagrangian of baryon interaction with the lepton field can be written in 



^In what follows we use the system of units h = c = 1 and the Boltzmann constant kB = 1- 




Ib {sb) 
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(1) 



exp {{sb + Ub- fiB) IT) + 1 
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a point-like approach (summation over repeated Greek indexes is assumed) 

■'-weak J V / 

where Gp is the Fermi weak couphng constant, and the Cabibbo factor C = cos9c = 0. 973 
for change of strangeness AS* = and C = sin 6'c for A^* = 1. For the direct Urea process 
involving baryons Bi and B2 

Bi^B2 + l + ui, (3) 
the lepton and baryon weak charged currents are, respectively: 

r = Ua'^ (1 - 75) Z/;, Ja = 1p27a{Cv - CaTs)^'!- (4) 

Here ipi and ip2 stand for the initial and final baryon fields; Cy and Ca are the corresponding 
vector and axial- vector coupling constants, respectively. In what follows we consider massless 
neutrinos of energy and momentum qi = {ui, ki) with Ui = ki. The energy and momentum 



of the final lepton of mass mi is denoted as q2 = {002, ^2) with 002 = y k2 + mf. 

We consider the total energy which is emitted into neutrino and antineutrino per unit 
volume and time. Within /3-equilibrium, the inverse reaction B2 + I Bi + ui corresponding 
to the capture of a lepton /, gives the same emissivity as the decay but in neutrinos. 
Thus, the total energy loss Q for the Urea processes is twice more than that caused by 
/?-decay (|^). Taking this into account by Fermi's "golden" rule we have 

^_^ Gj,C^ f d'k2d^kld^P2d'pi 2, , f .w. .^ 

X (27r)^ 6 {ei + Ui - 62 - U2 - uji - U2) S (pi - p2 - ki - kg) , (5) 

where the square of the matrix element of the reaction (j^) summed over spins of initial and 
final particles has the following form 

\Mf,f = 64 [(Cl - C^y) M*M; (giga) + {Ca - Cyf (giP^) (92^1) 

+ (C^ + Cy)'(giPi)(g2P2)l (6) 



with Pi = (£i,pi) and P2 = (e2,P2)- Antineutrinos are assumed to be freely escaping. The 
distribution function of initial baryons Bi as well as blocking of final states of the baryon 
B2 and the lepton / are taken into account by the Pauli blocking-factor /i (1 — /2) (1 — /O 
with 

fi M = J (7) 

exp {uj2 - M /T + 1 

being the Fermi-Dirac distribution function of leptons with the chemical potential fii. By 
neglecting the chemical potential of escaping neutrinos, we can write the condition of chem- 
ical equilibrium as fii = ^1—^2- We assume degenerate matter. Then by the use of the 
energy conservation equation ei + Ui = 62 + U2 + uj2 + ^^i, and taking the total energy of 
the final lepton and antineutrino as uj2 + uji = fii + uj' we obtain the identity 

/l (£1) (1 - f2 (62)) (1 - fl M) 

^ h (£1) (1 - h (^1 - ^0) (1 - fl if^i + ^' - ^1)) , (8) 

where uj' ^ T because the energy exchange in the direct Urea reaction goes naturally on the 
temperature scale ~ T. Due to the strong degeneracy of the medium, the main contribution 
to the integral @ comes from narrow regions of momentum space near the corresponding 
Fermi surfaces. Since the neutrino energy is uji ~ T, and the neutrino momentum fci ~ T is 
much smaller than the momenta of other particles, we can neglect the neutrino contributions 
in the energy-momentum conserving delta-functions. Then 

5 {si + Ui - 82 ~ U2 - uji - 002) S (Pi - P2 - ki - k2) 

-6 (cose - -V (pi -pI + kl)] 6 (pi - P2 - k2) , (9) 



Pik2 \ 2pik2 ^ ' ) 

where the Q is the angle between the momentum pi of the initial baryon and the momentum 
k2 of the final lepton. Now integrations over p2 momenta and all solid angles can be easily 
done. As mentioned above, we can replace the momenta of particles by the corresponding 
Fermi momenta, i.e. pi = "ppx^ p2 = Pf2, ^2 = Pfi in all smooth functions of energy and 
momentum under the integral. Then the remaining integration is reduced to the following 



duJiujldJ dei /i (^i) (1 - /i (ei - uj')) {I - fi{^ii + uj' - uji)) 



rfw' -, / duji -— = TT^T^ (10) 

expw'/T-lio 1 + exp (wi - w') /T 5040 ^ ^ 

Finally, to the lowest order in T//i; we obtain the following neutrino emissivityQ: 

Q = Y^^fC'T' [CvCa ((e^i + eF^)pli - (e^i - ep^) {pi, -pi,)) 
+2ClfiiM*M; + (C^ + d) [fil {2epieF2 - M*M;) + sp^pl, 
-\i.<^Fi + eF2)(pli-pl2 + Pli)) Q {Pfi + Pf2 - Pfi) (11) 



with G (x) = 1 if X > and zero otherwise. When the baryon and lepton momenta are at 
their individual Fermi surfaces, the 6- function (^ contributes only if pf2 +Pfi > Pfi- This 
"triangle" condition required by the step- function in Eq. (|Tl|) defines the above mantioned 
threshold for the direct Urea reactions. 

Eq.([TT|) can be obtained also by the use of the medium response functions. By the 
fluctuation-dissipation theorem, the total emissivity can be written as: 

GlC^ f il-fiiuj2))2lm[wfL^p\ d^k, d^h 



2 J exp(go-/ii + /i2)/T-l (27r)3 (27r)3 ^ 

where the integration goes over the phase volume of the final lepton and antineutrino of 
total energy go = ^^i + ^2 and total momentum q = k^^ + k2. The final-state blocking of the 
outgoing lepton is taken into account by the Pauli blocking-factor (1 — // (co'2)). The factor 
[exp (go — fJ'i + (^2) /T — l]^^ in Eq. (|T2]) arises due to the fact that the baryons Bi and 



B2 are in thermal equilibrium at temperature T and in chemical equilibrium with chemical 
potentials /xi and /X2, respectively. L"^ is defined here by 



a/3 



{kl4 + k'^k^, - {k, ■ k,) - le'-"''' (k,)^ {k2)s) , (13) 



^ We have corrected two misprints appearing in the published version of our paper (see Eq. (11) 
in ref. M). 
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and W^^ is the retarded weak-polarization tensor of the medium. The imaginary part of 
the retarded polarization is related to that of the causal (or time ordered) polarization Wap 
as follows 

Im Wf = tanh Im W^^. (14) 

Since the baryon weak-current includes vector and axial-vector contributions, to the lowest 
order in weak interactions, the polarization tensor consists on those diagrams which have 
ends at the weak vertex (6*1/7" ~ 6*^7^75) • Thus, we can take the weak polarization tensor 
as the sum of vector-vector, axial-axial, and mixed axial-vector pieces: 

= C^W^ + Cllif - 2CvCaIIva- (15) 
The vector-vector polarization tensor is of the following general form 



= ~^J '^-^ Tr [ G, (P° - t/i + go, P + q) I'^G, (P° - U,, p) 7^^ 



(16) 



were = (go, q) is the total four- momentum transfer, and the baryon propagators Gi and 
G2 correspond to different baryons Bi and B2. In the following we consider a not too large 
momentum transfer and not too high temperatures, permitting us to neglect antibaryons in 
the system. Then the in-medium baryon propagator can be written as: 



Gb{P) = {iPb + M*^) 



(17) 



The first term in Eq. (|T7|) describes propagation of free baryons, and the second term includes 
the Pauli principle restrictions^. Thus, the polarization tensors is the sum of a Feynman 
piece and an explicitly density-dependent piece. The energy-momentum transferred in the 



^Notice that this form of the baryon propagator takes into account not only the Pauli principle, 
but also partially strong interactions, which are included by means of the effective mass of the 
baryon and the self-consistent nuclear field. 
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considered processes are not too large (gg — — mf, with mi being the lepton mass), 
therefore we neglect the vacuum contribution and consider only the density-dependent part 
of polarization functions that are related to virtual excitations of a particle-hole type, with 

By a simple change of the integration variable Po = Pq — U2, Eq. ( [TE| ) can be reduced to 
the standard form 



n 



a/3 



/ 75lTr[Gi(p + g)7°G2(p)7^], (18) 



(2vr)^ 

with = {po, p) but with the effective momentum transfer given by 

0'^ = (go,q), % = qo-Ui + U2. (19) 

Obviously, the polarization functions depend explicitly on go instead of the total energy 
transfer go- The only known generalization (see refs. |^, @]) of the polarization functions to 
the case of different baryons in the polarization loop have been made by the ansatz that the 
general form of the tensor ([T8|) should be a sum of longitudinal and transverse components. 
Such an approach (which is valid in vacuum under the assumption of isotopic invariance of 
nucleons) is invalid for a system of interacting baryons, because in a medium the isovector 
current, caused by conversion of the baryon Bi ^ B2, is not conserved. 

Let's examine the nonconservation of the baryon current by contracting together the 
vector-vector polarization tensor and the effective four-momentum transfer 

Q„n-^ = _z I ^ Tr [G,{p + Q) [^Q] G^^ph^] . (20) 

The baryon propagator (P^) obeys the equation {'yp — M"^) Gb{p) = 1 with 1 being the 
identity matrix. By taking this into account and using the following identity 

go7o - q7=[(Po+go) 7o - (P + q) 7 - ^il - [Polo - P7 - M*] 

+ (M* - M*) (21) 



we obtain from Eq. (H 



+ (M*-M*) n^, 



where 11"^ is the mixed vector-scalar polarization, given by 



(22) 



— ? 



Tr 



Gi(p + Q) G2(p)7^ 



(23) 



We can make the change p + Q — > p in the second integral of Eg. (p2D . Then, using the 
explicit form ([T7| ) of the baryon propagators, we obtain 



Q^U"" = 6po {n2 -n,) + (M* - M*) U^, 



(24) 



where ni and n2 are the number densities of baryons Bi and B2, respectively, and Spo = 1 
if (3 = and zero otherwise. Eq. (^) can be equivalently written as 



0/3 



V (ni-n2) + (M*-M*) H^. 



(25) 



Thus, in a medium, the baryon vector current is conserved only in the case of symmetric 
nuclear matter [rii = 77,2) and isotopic invariance of the baryons (M* = Mg), which together 
with Ui = U2 provide g°n°^-g^m^ = 0. If at least one of the above conditions is not fulfilled, 
then Qall"^ 7^ 0. In particular, this means that the vector-vector polarization tensor for 
charged currents can not be written in terms of longitudinal and transverse components, 
as assumed in refs. ( [Q, 0). In a frame of reference where Q = (go,0,0,g), the vector- 
vector polarization tensor (p!8| ) has four independent components 11°°, 11°^ = 11^°, 11'^^, and 

nil ^ ^22 

The axial-axial tensor 



^'a (Qo, q) = 
as well as the vector-axial 



d'^p 



Tr 



^1(1) + (5)7^75^2(^)7^5 



(26) 



Gi{p + Q)r ^2(^)7^5 



(27) 



are of the standard form: 



(28) 



(29) 



Our calculation of the integrals yields the following imaginary components of the retarded 
polarization tensors: 



Here 



imn$° 



1 



27iq 



(30) 



Imn^2 = imng' 



+ J (t - Mf + Mf) /o 



27rg2 



go 



(31) 



imnf = 



1 



1 



g2 ((M* - M*f - ql) + (g^ _ Mf + M*^) j 



(32) 



imn]j^ = Imnf = 



t/2 + go - + Mf ) /i 



(33) 



im n 



R 



27iq 



M*M*Io, 



imn^^ = 



87rg2 



2th + {t- Mf + ) 



(34) 



. _ ~2 2 

t = qo-q , 



(35) 



and the functions /„ are defined as 
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4 = tanh j £ de2 £^ [(1 - /2 (^2)) (^2 + g^) 

+/2(e2)(l-/2(^2 + g;))]. (36) 

with 

^0 = ?o -/^i +/^2- (37) 

The lower cut-off in Eq.(PBD arises due to Icinematical restrictions and reads 

^o = -|(t-Mi*2 + Mf) 

- (Ml + Aqf) [t - (M| - MlY). (38) 

This energy is real- valued if either t > (M^ + M^f or t < (Mj* - Ml)^ The first case 
corresponds to creation of a real baryon pair in the medium and requires inclusion of the 
vacuum (Feynman) piece of polarizations for a correct calculation. As mentioned above, we 
do not consider so large values of t. Thus, the imaginary part of the polarization functions 
does not vanish for t < (Mj* — M^)^. In terms of go and q, this can be written as 

(go - f/i + f/2)' < g' + {Ml - M;f . (39) 

Contrarily to the statement made in refs. ( [0) that the imaginary part of polarizations 
vanishes when the momentum transfer is time-like, our calculation shows some domain of 
go > g where the imaginary part does not vanish even if we assume = M| (in fact, one 
can neglect the neutron-proton mass difference in nuclear matter). This result is obviously 
correct, otherwise, the neutron decay process n ^ p + + u , which requires go " g^ > 
for final leptons, would not be possible in neutron stars. 

The nuclear matter in the core of neutron stars becomes transparent for escaping neu- 
trinos only at relatively low temperatures. Therefore we focus on the low-temperature limit 
(yUfi/T ^ 1) and consider degenerate baryons and leptons under thermal and /3-equilibrium. 
In this case the kinematical condition (|38| ) restricts possible values of the momentum trans- 
fer. Obviously, in the degenerate case the cut-off ( |58| ) of the integral ( |5BD should be lower 
than the baryon Fermi energy 
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eF2 >-f {t- Mf + Mf ) 



2t 



- (M| + Mlf) {t - (M| - Mtf). 



(40) 



As mentioned above, the energy of the final lepton is close to its individual Fermi energy 
fii, while the condition of chemical equilibrium is /iz = /ii — /i2, where the baryon chemical 
potentials can be approximated by their individual Fermi energies. Thus, to the lowest 
order in T/ fii we can take go = ^fi — £f2, and t = {efi — SF2f — <f in all smooth functions. 
Inserting this into Eq.(HO) we obtain the following domain of the momentum transfer 



Pfi -Pf2< q< Pfi + PF2, 



(41) 



where the imaginary part of the polarization functions does not vanish. 
To the lowest accuracy in T/fii, we obtain from (RO): 



Jo ^ %, h ~ goeF2, h = q'oep^, 

where gQ = go "~ + is of the order of T. 

By taking also go = — in all smooth functions of Eqs. (|3 
(PI) we obtain 



with 



Svrg 



imn 



VA % % 



R 



inq 



im m 



57rg 



(42) 



in the interval 



(43) 



00 



[eFi + eF2Y-q'-iM*-M; 



(44) 



$03 = $30 



(^Fi - eF2) { {sfi + eF2f - g^ 



{eFi + eF2) {Mf-M;' 



(45) 



$33 = ^ 



[sfi - eF2f ([sfi + eF2f - 



If 



-2 (4i-42) (m*2-m^ 

+g2 (M* - M;f + (m*2 - Mf 



(46) 
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$11 = $22 = ^ 



1 

2g2 



(47) 



(48) 



= ^ ((e^Fi - £^F2) (pfi - - (^Fi + £f2) q^) ■ (49) 

Contracting together the tensors, as indicated in Eq. ([T^) , and calculating the imaginary 
part to the lowest order in T//i; we obtain 



Im 



=>33 



(50) 



+2go$" + - 2go$°^ + - 2g<l>") cos0 

-2C\ (go - q cos e) $^ - ACvCau^i [q - qo cos 6) <!>va\ ^(go, g), 

where 6 is the angle between the ki and q momenta. The step-function 

^(go, ?) = ©('? + Pf2 - Pfi) (pfi + Pf2 - q) , 
with (x) = 1 if z > and zero otherwise, comes from the kinematical restrictions (^l]) and 



(p9|). By inserting expression (^) into Eq. and taking go — /^z, ? — PFh we arrive to 

the integral over cigg and the neutrino phase space. We omit this straightforward calculation. 



which results in the expression given by Eq. (pUf). We stress here the fact that agreement 
between these two methods is achieved only when we take into account non-conservation of 
the baryon vector current. 

In a non-relativistic limit {Pfb ^ Af^, -C M%), the leading terms of Eq. ( pA]) give 
the following expression 

4577r 



Ql = Y^G'fC {Cy + 3^1) M^M^i^iT'e ippi + PF2 - PFi) 



(51) 



which coincides with the known non-relativistic result of Lattimer et al. M. As in that case. 



our formula ([Tl|) exhibits a threshold dependence on the proton concentration. The "trian- 
gle" condition Pfi+Pf2 > Pfi, required by the step-function, is necessary for conservation of 
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the total momentum in the reaction. We pay attention, however, that this condition should 
be supplemented by the one of chemical equilibrium. As mentioned in introduction, in the 
relativistic regime, the chemical equilibrium is possible only due to the fact that strong 
interactions create a gap between the proton and neutron energy spectrums, that is much 
larger than the mass difference of participating baryons. 

The relative efficiency of the direct Urea processes involving different kinds of baryons 
depends essentially on the composition of the /5-stable nuclear matter in the core of a neutron 
star. In order to quantify the relativistic effects in the direct Urea processes, we consider 
a simplified model for degenerate nuclear matter of the standard composition consisting on 
neutrons, protons and electrons under beta-equilibrium. We use a Walecka-type self- 
consistent relativistic model of nuclear matter, by assuming that the interactions among 
nucleons are mediated by exchange of a, uo, and p mesons. In this model the baryon effective 
mass M* = M — g^sO'o as well as the potential energies Ub = gujE^o + gpstsBbo are calculated 
in a self-consistent way. Here ujq, Bq, and o"o are, respectively, the mean-field values of the 
u, p, and (T-mesons; g^^B, QpBi and g^B are the strong interaction coupling constants, and 
t^B is the third component of isospin for the baryons. Parameters of the model are chosen, 
as in ref. [0, to reproduce the nuclear matter equilibrium density, the binding energy per 
nucleon, the symmetry energy, the compression modulus, and the nucleon Dirac effective 
mass at saturation density = 0.16 fm~^ . We have plotted in Fig. 1 the emissivity of the 
Urea process, as given by our formula Eq. (|ll]) in comparison with the emissivity ( pi]) given 
by Lattimer et al. Both magnitudes are normalized with respect to the emissivity given 
by our formula at threshold density. We can observe that relativistic effects dramatically 
modify the emissivity. The non-relativistic approximation approaches our result only at 
densities much smaller than the threshold density, indicated by an arrow on the density 
axis. Due to the decrease of the nucleon effective mass, the formula derived by Lattimer 
et al. predicts a decreasing of the emissivity as density increases above the threshold. In 
contrast, our formula shows a substantial increasing of the neutrino energy losses as we go 
to larger densities. Of course, the exact value depends on the underlying nuclear matter 
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model. 

Our Eq. (|lT]) is obtained in the mean field approximation and does not take into account 
correlations among the baryons, which normally have a tendency to suppress neutrino energy 
losses. These calculations are out of the scope of this letter, and we leave these unwieldy 
calculations for a future work. At some densities neutrons may be superfiuid and/or protons 
superconducting. The direct Urea rate is then reduced. This can be taken into account in 
the standard way (See e. g. ^j). 
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Fig 1. Neutrino emissivity due to the Urea processes. We plot, as a dashed hue below the 
threshold density nth = l.TQrio, and as a solid line above rith our result, given by Eq.(|ll|). 
We also plot, for the sake of comparison, the emissivity given by Lattimer et al, as a dotted 
line below rith, and as a dashed-dotted line above nth- Both magnitueds are normalized 
with respect to the emissivity given by our formula at threshold density. We labeled the 
threshold density by an arrow on the abscisa axis, which shows the number density in units 
of saturation density no = 0.16/m~^. 
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